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RESEARCH  OBJECTIVES 

The  research  objectives  of  the  Princeton  Catalytic  Combustion  Program  for 
the  period  8/81  to  7/83  have  been  both  the  formulation  and  testing  of  a  tran¬ 
sient,  two-dimensional  catalytic  combustion  model  capable  of  predicting  the 
light-off  transient  of  catalytic  combustors  as  well  as  the  testing  and  eva¬ 
luation  of  a  platinum-doped  perovskite  catalyst  in  terms  of  its  high  temperature 
durability  and  low  temperature  light-off  characteristics. 


STATUS  OF  RESEARCH 


The  potential  of  catalytic  combustion  for  application  to  stationary  gas 
turbines  [1-7],  aircraft  gas  turbines  [8-13],  highway  vehicle  turbines  [14-18], 
and  boilers  [6]  has  now  been  widely  recognized.  The  utility  of  catalytic  com¬ 
bustors  rests  upon  their  ability  to  operate  efficiently  over  a  much  wider  range 
of  fuel  air  ratios  than  are  imposed  by  the  flammability  limits  of  conventional 
gas  burners.  The  catalytic  combustor  also  has  the  advantages  of  low  NOx  levels, 
uniform  exit  temperature  profile,  and  combustion  stability  over  broader  ranges 
of  equivalence  ratio. 

It  could  be  argued  convincingly  that  catalytic  combustion  is  ready  for 
extensive  commercial  use,  however,  areas  in  which  further  progress  could  be  made 
can  also  be  identified.  Justifiable  objectives  for  current  research  include 
lowering  the  ignition  temperature  and  at  the  same  time  increasing  durability  at 
high  temperature,  reducing  fatigue  due  to  thermal  stresses  during  transient 
operations,  and  developing  design  optimization  tools  to  isolate  promising  ranges 
of  the  innumerable  parameters  (see  Table  1).  The  first  is  a  problem  of  physics 
of  solids  but  the  second  and  third  are  questions  that  can  be  addressed  within 
the  thermal  sciences. 

Two  Dimensional,  Transient  Catalytic  Combustion  Model* 

When  honeycomb  monoliths  are  used,  as  is  often  the  case,  the  flow  in  the 
channels  is  generally  laminar  and  convection  and  diffusion  of  mass,  momentum  and 
energy  in  the  gas  and  thermal  diffusion  in  the  solid  can  be  characterized  quite 


•This  work  has  been  submitted  for  publication  to  The  Journal  of  Heat  and  Mass 
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accurately.  Elementary  homogeneous  and  heoterogeneous  reactions  cannot,  but 
global  mechanisms  of  engineering  usefulness  can  often  be  determined.  The  con¬ 
servation  equations  for  the  gas  and  the  solid  represent  all  these  processes  and 
allow  one  to  study  their  interactions  as  well.  In  their  simplest,  one  dimen¬ 
sional  form  such  equations  have  already  been  extensively  and  fruitfully  used 
[19,20].  In  particular,  within  the  context  of  transient  operations,  notable  is 
the  work  of  T'ien  [21]  that  revealed  many  of  their  essential  features. 

One-dimensional  models  have  the  considerable  advantage  of  simplicity  but 
are  inherently  limited  in  their  resolution  capabilities.  They  allow  one  to  com¬ 
pute  axial  gradients  through  the  use  of  global  transfer  functions  for  the 
controlling  radial  gradients.  These  functions  can  be  very  accurate  for  fully 
developed  steady  flows  but  require  significant  adjustments  for  transient 
entrance  flows.  For  design  optimization  studies,  two-dimensional  models  offer 
greater  accuracy,  since  they  allow  one  to  resolve  both  axial  and  radial  gra¬ 
dients,  at  greater  computational  costs.  Through  them  one  can  compute  also  the 
thermal  radial  and  axial  gradients  in  the  solid  that  are  relevant  to  the  problem 
of  thermal  fatigue.  Our  earlier  two  dimensional  model  [22-23]  was  only  for  the 
steady  state  of  the  gas  within  a  channel  of  the  monolith.  It  used  measured  wall 
temperatures  instead  of  solving  for  the  thermal  gradients  in  the  solid.  Thus, 
it  was  mostly  a  tool  for  the  detailed  interpretation  of  measured  quantities.  In 
the  work  reported  here  the  complete  coupled  two-dimensional  transient  equations 
for  both  solid  and  gas  are  solved,  however,  to  limit  computation  costs,  the 
reported  solutions  were  obtained  using  the  quasi-steady  gas  assumption  already 
employed  by  T'ien  in  his  one-dimensional  model. 

The  reported  work  concerns  the  oxidation  of  mixtures  of  CO  and  moist  air  on 
platinum.  The  oxidation  of  CO  is  an  important  step  in  the  combustion  of  hydro¬ 
carbons.  Platinum  was  chosen  because  it  is  one  of  the  best  characterized  cata- 
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lyst  even  though  it  is  not  the  best  one  for  catalytic  combustion  due  to  its 
limited  high  temperature  durability*  However,  its  low  temperature  ignition  pro¬ 
perties  are  excellent  and  it  remains  a  good  candidate  for  certain  applications. 

In  the  following  sections,  the  time  scales  associated  with  the  problem  are 
considered  first.  Then  the  equations  of  the  model  are  presented  and  a  detailed 
discussion  is  given  of  one  ignition  transient.  Finally  steady  state  results 
computed  through  the  transient  are  compared  with  measured  ones  (no  transient 
measurements  are  available).  Details  can  be  found  in  [24]. 

Characteristic  time  scales 

For  a  well-insulated,  monolithic  catalytic  combustor  with  uniformly  spaced 
channels,  the  study  of  the  entire  reactor  can  be  reduced  to  the  one  of  a  single 
channel.  Such  a  channel  is  shown  in  Fig.  1  and  is  considered  here.  The  details 
of  the  specific  catalyst  considered  are  in  Table  1.  The  characteristic  time 
scales  of  this  problem  are  listed  in  Table  2,  following  T'ien  [21]. 

In  computing  the  estimates  of  Table  2,  averaged  quantities  were  used.  This 
precludes  axial  variations.  Actually  the  time  scales  vary  along  the  length  of 
the  channel.  For  example,  in  the  entrance  zone,  the  boundary  layer  growth  is 
very  fast  and  consequently  the  rates  of  transport  vary  rapidly,  as  well.  Also, 
local  sources  and  sinks  of  heat  reduce  local  length  scales.  Thus  the  times  of 
Table  2  are  very  rough  estimates.  Indeed,  the  substrate  heat  up  time  will  ulti¬ 
mately  be  computed  to  be  between  1  and  10  s  instead  of  60s.  Nonetheless  the 
estimates  suggest  some  important  simplications  whose  validity  can  always  be 
checked  aposteriori. 

The  characteristic  times  of  the  gas  are  shorter  by  about  one  order  of 
magnitude  than  those  of  the  solid  so  that  the  quasi -steady  gas  phase  assumption 
may  be  justified  in  our  case.  However,  the  difference  between  gas  and  solid 
times  is  not  overwhelmingly  large  and  is  sensitive  to  geometrical  parameters. 


TABLE  1.  CATALYST  PROPERTIES 


SUBSTRATE:  Cordierite,  American  Lava  Corp. ,  AlSiMag  795,  split  cell 

length  0.076Q£  0.003  m 

wall  thickness  0.25  x  10-3  m 

open  area  64% 

channels  per  unit  area  0.35  x  106  m-2 

open  area  per  channel  1.9  x  10-6  m2 

ratio  surface  area  to  total  volume  2100  o2/m3 

ratio  surface  area  to  gas  volume  3260  m2/m3 

channel  hydraulic  diameter  0.0014  m 

bulk  density  610  kg/m3 

solid  density  1 700  kg/m3 

safe  operating  temperature  1473  K 

specific  heat  800  J/kg* K 

coefficient  of  thermal  expansion  (linear,  294-1033  K)  3.8  x  10-6  K-1 
thermal  conductivity  of  solid  (572  K)  1.4  J/m*  s»  K 

approximate  channel  cross-section  is  a  trapezoid  with  base  lengths  .0011  m 
and  .0023  m,  and  height  .0012  m 
WASHCOAT:  y -alumina 
loading  115-125  kg/m3 

surface  area  (BET) (29. 2-33. 0)  x  106  m2/m3 
CATALYST:  platinum 
loading  4.2  kg/m3 

surface  area  (CO  chemisorption)  6  x  1 04  m2/m3 

They -alumina  and  platinum  were  applied  by  Matthey  Bishop,  Inc.,  Malvern,  PA. 
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TABLE  2 


ESTIMATE  OF  TIME  SCALES  OF  A  CATALYTIC  CHANNEL  IN  THE  MONOLITH  OF  TABLE  I 


1 

TIME 

DEFINITION 

MAGNITUDE 

i 

i. 

Gas  residence  time 
in  reactor 

l/°in 

(1.5  -  7.5)  x  10~3  s. 

!■ 

I* 

2. 

Heat,  mass  &  momentum 

transfer  time  between 

gas  and  catalyst 

surface 

Nu  a 

g  g 

1°"3  s 

D  t 

|« 

3. 

Homogeneous  {k  exp(- 

reaction  time 

V*T)I02,I820|,/21 

Variable  depending  on 
local  fuel  concentra- 

tion  and  temperature 

4. 

k 

Heterogeneous  { 2  — 

reaction  time 

eXp(-Es/RT)^"1 

Variable  depending  on 
local  fuel  concentra- 

tion  and  temperature 

IQ 

5. 

Time  for  temp¬ 
erature  wave  to 

(s/2 ) 2 
a 

s 

10-2  s. 

reach  the  substrate 

half  thickness 

6. 

Time  for  temperature 
wave  to  reach  the 

a 

g 

3  x  10-3  s. 

centerline  of  the 

i( 

gas  channel 

7. 

Time  for  temperature 
wave  to  reach  the  end 

L^ 

a 

s 

5  x  1 03  s . 

of  the  substrate 

‘ 

8. 

Substrate  heat-up 
time 

i  f>  cp_ 

1  s  s  R 

Nu  k  s 

g  g 

60  s. 

10-3  -  10-1  S. 


( 


9.  Input  characteristic 
time 
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Thus  the  quasi-steady  gas  approximation  cannot  be  assumed  valid  for  all  cataly¬ 
tic  combustors  of  practical  interest.  Similarly,  in  the  solid,  the  time  for 
radial  thermal  equilibrium  is  much  shorter  than  that  for  axial  equilibrium  and 
one  would  be  justified  in  treating  the  solid  as  one  dimensional,  however  the 
important  thermal  stresses  would  not  be  properly  resolved.  Our  model  does  not 
require  either  one  of  the  above  two  limitations,  nonetheless,  the  reported 
results  were  obtained  with  the  quasi-steady  gas  phase  approximation  to  reduce 
computation  costs.  Radial  gradients  in  the  solid  were,  however,  resolved. 

Notice  that  associated  with  heterogeneous  fuel  conversion  are  two  time 
scales:  that  for  the  surface  reaction  and  that  for  the  transfer  of  mass  between 
the  gas  and  the  surface.  When  the  former  is  much  smaller  than  the  latter  fuel 
conversion  is  mass  transfer  limited.  This  is  never  the  case  at  the  beginning  of 
the  transient  but  is  often  the  case  in  steady  state. 

Equations 

The  gas  phase  equations  for  the  flow  inside  a  channel  schematized  in  Fig. 

1,  are  the  complete,  two-dimensional,  conservation  equations  for  reactive,  lami¬ 
nar  flow  [25].  For  the  experimental  conditions  investigated,  the  Reynolds 
number  of  the  flow  is  below  the  critical  one  (~  2000)  for  transition  to  tur¬ 
bulence. 

The  general  conservation  equation  for  the  gas  is  written  as 

f£l*i  ♦!_  (prv*)-|j(rr|t)  -Ijfrrlt,, 


1 


« 


where  $ 


r 


S 


(1 ) 


=  u,  v,  T,  y^,  etc. 

=  appropriate  transport  property,  such  as  p,  K,  pD,^,  etc. 

=  corresponding  source  term,  such  as  -3p/3x  or  the  rate 
of  species  formation/destruction 


3<f>  /3 1 


unsteady  or  inertial  term 


3p  /3x 
3p  /3  r 

3  (rr  3f/5x)/3x  =  axial 

}  diffusive  terms 

3  (rr  2<f>/3r)/9r  =  radial 

For  the  substrate  only  the  energy  conservation  equation  is  necessary.  It 
has  the  same  form  as  the  general  conservation  Eq.  1  except  for  exclusion  of  the 
convective  and  source  terms: 


=  axial 


}  convective  terms 


=  radial 


3  T  1  .3 

pwcp  3t  "  r  3 x 
w 


(r 


K  51  ) 
v  3x  1 


3  .  v  3T  . 

57  (r  KW37  )] 


(2) 


AH  the  model  equations  are  given  explicitly  in  Table  3.  The  laminar  tran¬ 
sient  coefficients  (p ,  K,  Dik )  and  the  specific  heats  varied  from  species  to 
species  and  with  temperature  and  were  obtained  from  Refs.  26  and  27  respec¬ 
tively. 

The  approach  utilized  to  solve  the  equations  may  be  called  a  "continuum" 
approach.  The  solution  algorithm  makes  no  distinction  between  the  gas  phase  and 
the  substrate.  The  material  properties  change  appropriately  and  the  convective 
terms  are  set  to  zero  for  the  substrate.  Temperature  matching  at  the  interface 
is  specified  as  a  boundary  condition.  With  this  approach,  the  gas  and  the 
substrate  energy  equations  are  solved  simultaneously  and  temperature  matching  at 


the  interface  is  automatic 


I 

Boundary  and  initial  conditions 

a)  The  boundary  conditions  used  in  the  model  are: 

(i )  At  the  inlet  (x=0)  of  the  gas  phase,  all  physical  variables  are  spe 
cified  including  the  gas  pressure  and  velocity.  The  fuel  concentration 
increases  with  time  until  it  reaches  its  final  steady  state  value  according  to 


WT  > 


CO,  F 


[1 


(3) 


and  for  YcO*T  ^  >  0.99  YC0/p  »  Yco^T  ^  ~  *CO,F 


where  t 
t 

t* 


YC0(T  ) 
yCO,F 


P 


non-dimensional  time  =  t/t* 
time 

characteristic  variable  for  fuel  ramp 

~  10~3  -  10-1  sec. 

fuel  concentration  at  inlet 

fuel  concentration  at  inlet  at  the  end  of  the 
fuel  ramp. 


(ii)  At  the  outlet  (x=L)  of  the  gas  phase,  all  physical  variables  4>  have 
their  gradients  set  to  zero,  i.e.  3$ /) x  =  0.  This  linear  condition  simplifies 
the  solution  and  is  also  observed  in  practice,  provided  the  L/D  ratio  of  the 
channel  is  (  ~  100). 


(iii)  For  the  substrate,  the  only  variable  of  interest  is  the  temperature 
T.  The  axial  temperature  gradients  must  satisfy  the  radiative  heat  exchanges 
from  the  end  faces  of  the  catalyst  to  the  upstream  and  downstream  parts  of  the 


test  section: 


conductivity  of  the  substrate 


where  K 
w 


e  =  emissivity  of  the  substrate 


o  *  Stef an-Boltzmann  constant 


T  ,  T  =  temperature  of  the  substrate  face  at  x  *  0  and 

f,o  f,L  x  =  L  respectively 


tTS,IN'  tTS,OUT  =  temperature  of  the  test-section  at  the  inlet 

and  the  outlet  respectively  =  600  K 

(iv)  At  the  centerline  of  the  channel  (r=0),  the  boundary  conditions 
express  radial  symmetry,  i.e.  3<|)/3r  =  0. 

(v)  At  the  half-thickness  of  the  substrate  (r=R+s/2),  small  radial  losses 
are  present  in  the  experiment  in  the  magnitude  of  10  K/channel  [23].  To  account 
approximately  for  them,  the  following  temperature  gradient  is  specified: 


dT 

dr  r  =  R  +  ^ 


tnjwmi  "  tnjw 
— 


(5) 


where  TNJW 

tNJWMI 


Lx 


temperature  at  the  substrate  half -thickness 
temperature  at  the  mesh  point  next  to  the 
substrate  half -thickness 
local  radial  cell  length. 
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Assume  (TNJW  -  TiN)  =  C(TNJWMI  -  TiN) 


dr  I  r  =  R  +  -^ 


(1-C) 


NJWMI 

Ar 


Thus,  the  simulation  of  radial  losses  is  reduced  to  the  specification  of 
the  single  constant  C.  The  value  of  C  was  0.9985. 

(vi )  At  the  gas-substrate  interface  (r=R),  the  no-slip  conditions  for  the 
gas  apply,  i.e.  u  =  v  =  0  and  the  temperature  satisfies  the  balance  between  the 
catalytic  heat  release  at  the  surface  and  conduction  into  the  substrate  and  the 


gas  phase 


'‘“gUV  -  ‘  ^catalytic 


where  Kg  is  the  conductivity  of  the  gas,  Ks  is  the  conductivity  of  the 
substrate,  and  ^catalytic  is  tfie  rate  heat  release  at  the  wall  per  unit 


The  chemical  species,  N2  and  H20  are  assumed  inert  and  their  radial  gra¬ 
dients  set  equal  to  zero.  When  the  reaction  at  the  wall  is  not  diffusion 
controlled,  a  kinetic  rate  is  imposed  for  CO. 


o  D  i  s  m 

p  CO,N23r  r=R  CO 


where  mCQ  is  the  CO  mass  consumption  rate  per  unit  area. 


When  the  process  gets  diffusion  limited,  a  switch  to  infinitely  fast  kine¬ 
tics  is  made  automatically  and  the  boundary  condition  for  CO  becomes  Vco  = 


Further  assuming  that  at  the  wall,  CO  is  oxidized  via  the  reaction 


CO  +  1/2  O2  ♦  co2,  the  boundary  conditions  for  O2  and  CO2  are 


3*. 
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b)  The  initial  conditions  correspond  to  the  situation  when  there  is  no 
fuel  in  the  system.  Non-reacting  air  at  the  initial  temperature  Tj^  flows 
through  the  channel  and  the  entire  system  is  in  equilibrium  at  this  temperature. 
Therefore 


T(x,r,0)  =  Tin 

Yi(x,r,0)  =0  i  =  CO,  C02  (9) 

Yi(x,r,0)  =  V.fIN  i  -  02,  N2,  H20. 

Heterogeneous  reaction 

A  vast  literature  [28-32]  exists  on  the  oxidation  of  CO  on  platinum.  But 
the  available  rates  apply  to  steady  state.  In  the  present  study,  one  deals  with 
a  transient  situation  in  which  the  fuel  concentration  at  the  surface  of  the 
catalyst  changes  rapidly.  A  qualitative  study  of  the  unsteady  kinetics  of  the 
heterogeneous  reaction  of  CO  and  oxygen  on  platinum  was  presented  by  Bykov  and 
Yablonskii  [33].  Conversion  rates  of  the  order  of  lO1^  moles/cm2 -s  were 
reported.  The  rates  rise  very  rapidly  with  fuel  concentration  to  these  values 


i 
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and  then  fall  off  sharply  to  the  steady  state  values.  There  is  a  rate  pulse 
which  triggers  the  catalytic  reaction.  However,  this  was  a  low  pressure  study 
(~10“6  Torr)  and  there  is  uncertainty  in  extrapolating  these  results  to  our 
higher  pressure  (~10+^  Torr)  situation. 

In  the  absence  of  any  concrete  data,  the  reaction  mechanism  and  rate 
egressions  chosen  are  simple  ones.  The  rate  constant  is  kept  sufficiently  high 
to  ensure  that  the  rate  of  temperature  rise  of  the  system  is  compatible  with 
that  observed  experimentally.  The  model  chosen  is  a  simple,  one-step  overall 
reaction  corresponding  to  the  reaction  betwen  adsorbed  fuel  and  adsorbed  oxygen: 


CO+2°2  +  C°2 


The  reaction  rate  is  taken  as: 


d  [CO] 
dt 


k  [CO]  exp (-E  /R  T ) 
s  s 


moles 

2 

cm  -s 


with  k  =  3  x  10  cm/s,  E  =  24.9  kcal/mole,  and 

s  s 

[CO]  is  the  concentration  of  CO  near  the  wall  in  moles/cm3. 

The  activation  energy  agrees  with  that  used  in  [32,34,35]  which  cite 
overall  rates. 

Homogeneous  kinetics 

The  model  chosen  for  the  homogeneous  oxidation  of  CO  by  O2  is  a  simple, 
one-step  overall  reaction: 


CO  +  i  °2  -  C02 


The  reaction  rate  is  taken  as  in  [36,23]; 


< 


d  [CO] 
dt 


k  [COHO  ]0,5[H  O]0,5  exp(-E  /R°T)  -??^-eS 
o  Z  /  O  J 

cm  -s 


Hi  > 


14  3 

where  k  =  1.3x10  cm  /mole-s,  E  =126  KJ/mole  =  30  kcal/mole, 

o  °  3 

and  [  ]  =  concentration  in  moles/cm 

Radiation  model 

The  test-section,  which  has  a  square  cross  section,  is  modeled  as  a 
cylinder.  Radiation  is  thus  simplified  to  the  calculation  of  radiative  losses 
from  the  front  and  back  surfaces  of  the  solid  and  the  inside  of  a  cylinder  (the 
channel)  to  the  inside  of  another  cylinder  (the  modeled  test-section).  The 
angle  factors  for  this  geometry  can  be  easily  evaluated.  The  model  computes  the 
radiative  loss  from  each  computational  element  on  the  surface  of  the  catalyst. 
Similarly,  radiation  inter-exchange  between  different  elements  on  the  catalyst 
surface  is  also  evaluated  [24].  In  our  case  the  net  effect  of  radiation  is  heat 
losses  that  reduce  the  energy  available  for  conduction  from  the  interface  into 
the  gas  and  the  solid. 

Outline  of  the  solution  algorithm 

All  details  of  the  solution  technique  are  available  in  [24],  including 
the  listing  of  the  program.  Here  only  a  brief  outline  is  given. 

The  equations  are  solved  numerically  by  a  method  based  on  the  TEACH  code 
documented  by  Gosman  and  Ideriah  [37].  The  code  used  in  the  present  work  dif¬ 
fers  from  the  original  TEACH,  that  was  for  steady  incompressible  and  nonreactive 
flows,  by  the  introduction  of  species  and  energy  conservation  equations, 
variable  density  and  the  equation  of  state  which  links  pressure,  density  and 
mass  fractions  and  the  extension  to  handle  gas  and  solid  simultaneously  and 
through  their  transients. 

The  solution  algorithm  is  as  follows.  The  general  Eq.  1  is  integrated  over 


the  elemental  cell  volume,  so  that  only  convective  and  diffusive  fluxes,  and 


sources  appear  inside  the  computer  code.  The  partial  differential  equations  for 
the  conserved  quantities  are  thereby  transformed  into  linear  finite  difference 
equations  which  nay  be  solved  iteratively.  The  integration  domain  is  divided 
into  cells.  For  each  variable  a  new  line  of  values  at  each  longitudinal  loca¬ 
tion,  x,  is  computed.  Marching  down  the  axis  of  the  channel,  the  entire  channel 
is  swept.  After  the  entire  mesh  has  been  swept,  the  total  amount  of  conserved 
quantities  is  computed  and  compared  with  the  entering  fluxes.  The  process  is 
repeated  iteratively  if  the  difference  is  above  a  set  cut  off,  typically  5%. 

A  hybrid  differencing  scheme  is  used  that  reduces  to  upstream  differencing 
for  the  convection  terms  if  the  absolute  value  of  the  cell  Pellet  number  >  2  and 
to  central  differencing  if  <  2.  The  grid  is  finer  near  the  entrance  and  the 
solid-gas  interface  and  expands  radially  and  axially  from  those  two  locations. 
Standard  tests  of  increasing  numerical  resolution  and  of  comparing  numerical 
results  with  close  form  solutions  for  simpler  cases  were  performed  and  are 
discussed  in  Ref.  24. 

Conplete  transient  solutions  for  the  gas  and  the  solid  were  obtained  with 
an  implicit  time  scheme,  time  steps  of  the  order  of  10“^  s,  and  about  100  itera¬ 
tions  per  time  step.  But  the  computation  time  was  too  long  for  parametric  stu¬ 
dies  and  the  quasi-steady  gas  phase  approximation  was  made.  As  explained  in  the 
introduction,  this  approximation  is  not  always  valid  for  practical  combustors 
but  is  acceptable  for  the  conditions  of  interest  in  this  paper. 

RESULTS 

The  general  sequence  of  the  computed  transient  events  is  as  follows. 
Initially  there  is  no  CO  in  the  channel.  The  CO  concentration  at  the  entrance 
of  the  channel  is  gradually  increased  to  its  final  steady  state  value.  (This 
fuel  ramp  occurs  in  a  time  which  is  short  in  comparison  to  that  required  for  the 
entire  system  to  reach  its  steady  state  configuration. )  Thus  at  the  beginning 
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the  mixture  is  very  lean  and  the  temperature  low.  Conversion  starts  at  the  sur¬ 
face  and  its  rate  is  controlled  by  the  rate  of  the  heterogeneous  reaction.  Due 
to  the  surface  activity  a  boundary  layer  is  formed  that  remains  thin  near  the 
entrance  but  grows  thicker  along  the  length  of  the  channel.  Initially  most  of 
the  heat  is  released  near  the  entrance  and  goes  to  increase  the  temperature  of 
both  solid  and  the  gas.  The  energy  received  by  the  gas  near  the  entrance  is 
convected  downstream  and  diffused  back  to  the  solid,  thus  increasing  its  back 
temperature.  Conduction  along  the  solid  occurs  simultaneously.  As  the  tem¬ 
perature  of  the  system  increases,  homogeneous  reactions  start  to  contribute, 
surface  reactions  become  faster  than  diffusion,  and  radiation  losses  become 
significant.  The  time  to  reach  steady  state  is  controlled  by  the  heating  of  the 
back  end  of  the  solid.  This  sequence  of  events  is  in  general  agreement  with 
that  discussed  by  T'ien  [21]. 

A  detailed  description  of  one  transient  is  now  given,  followed  by  a  com¬ 
parison  with  measured  substrate  temperatures  and  exhaust  emissions  from  a  series 
of  steady  state  experiments. 

Transient  study 

This  section  contains  the  results  for  a  case  with  the  following  final  con¬ 
ditions  at  the  inlet: 


a ) 

Velocity 

: 

10  m/s 

b) 

Temperature 

• 

600  K 

c) 

Pressure 

: 

110  K  Pa 

d) 

Equivalence  ratio 

: 

0.  32 

e ) 

Adiabatic  flame  temperature 

: 

1580  K 

f ) 

Water  mole  fraction 

: 

5.4  x  10 

g) 

Mach  number 

* 

0.02 

h) 

Reynolds  number 

5 

290.4 

( 


It  is  one  of  the  cases  of  Ref.  23  and  its  behaviour  is  typical  for  our 
system.  Initially  the  temperature  is  600  K  everywhere  and  there  is  no  fuel 
anywhere.  The  transient  is  begun  by  letting  fuel  into  the  system  and  the  fuel 
time-ramp  is  shown  in  Fig.  2. 

Overall  results 

The  following  energy  balance  can  be  written  for  the  system: 

rate  of  rate  of  rate  of 

(chemical  energy +  (kinetic  energy +  (thermal  energy ) 

=  (rate  of  chemical  energy )  t  +  (rate  of  kinetic  energy )QUt 

+  (rate  of  thermal  energy )QUt  +  (rate  of  radiation  loss  from 

within  the  channel)  +  (rate  of  thermal  energy  into  the  substrate) 
The  following  definitions  are  introduced  (hi  +  h 2  +  h 3  +  h 4  +  h 5  =  1 ) : 


mat 
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(rate  of  thermal  energy )QUt  "  (rate  of  thermal  energy 
(rate  of  chemical  and  kinetic  energy 


(rate  of  radiation  loss  from  the  channel) 
(rate  of  chemical  and  kinetic  energy). 


(rate  of  thermal  energy  into  the  substrate) 
(rate  of  chemical  and  kinetic  energy )^ 


(rate  of  kinetic  energy) 


(rate  of  chemical  and  kinetic  energy )^ 


(rate  of  chemical  energy )QUt 
(rate  of  chemical  and  kinetic  energy )~ 


Fig.  3  shows  the  energy  distribution  during  the  transient  and  represents  an 
overall  picture  of  the  transient  processes.  Notice  first  that  the  imposed  fuel 
transient  at  the  inlet  lasts  less  than  0.1  s  and  the  system  reaches  steady  state 
at  about  10s.  Thus,  the  results  do  reflect  the  transient  characteristics  of  the 
catalyst  and  not  those  of  the  fuel  delivery  system. 

Initially,  all  the  energy  release  goes  into  heating  the  substrate.  As 
shown  by  the  behaviour  of  n  3 •  This  is  true  until  t  ■  0.1s  when  the  final  fuel 
concentration  gets  established  at  the  inlet.  The  initial  heat  release  occurs  on 
the  surface  of  the  catalyst.  The  process  of  equilibrating  all  the  substrate  but 
its  entrance  begins  towards  the  end  of  the  fuel  ramp.  This  can  also  be  seen  by 
the  stabilization  of  the  radiation  losses  n 2  at  a  value  fairly  close  to  the 
final  one  (Radiation  losses  occur  primarily  from  the  front  of  the  catalyst).  At 


the  same  time  the  energy  content  of  the  gas  phase,  n 1 ,  is  increasing.  This  has 
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FIGURE  3.  DISTRIBUTION  OF  INCOMING  ENERGY  DURING  THE  TRANSIENT 
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two  effects.  Homogeneous  reactions  start  playing  a  more  important  role  and  very 
little  unburnt  fuel  finds  its  way  to  the  exhaust.  This  is  shown  by  the  drop  in 
ri5«  The  increase  in  the  thermal  content  of  the  gases  also  manifests  itself  in 
the  form  of  higher  kinetic  energy  and  is  shown  by  the  increase  in  114. 

The  term  H3  is  significant  in  that  it  gives  the  transient  time  of  the 
system.  At  steady  state,  the  net  energy  flowing  into  the  substrate  is  equal  to 
that  radiated  by  the  front  and  back  surfaces  of  the  channel  and  lost  radilly  to 
adjacent  channels.  In  this  case,  it  amounts  to  about  20%  of  the  chemical  and 
kinetic  energy  flowing  in. 

Fig.  (4)  compares  the  overall  conversion  by  the  heterogeneous  and  homoge¬ 
neous  paths.  The  quantity  being  plotted  is: 


fJ  Us  dA 


The  numerator  is  the  overall  rate  of  destruction  of  CO  due  to  homogeneous 
reactions,  while  the  denominator  expresses  the  corresponding  catalytic  conver¬ 
sion  rate.  As  indicated  by  the  earlier  plots,  during  the  early  part  of  the 
transient,  fuel  conversion  ■>  s  entirely  due  to  the  catalyst.  This  is  an  expected 
result  since  the  light-off  temperature  for  the  catalytic  reaction  is  low.  This 
situation  remains  unchanged  till  almost  steady  state  when  the  gas  temperature  is 
high  enough  for  the  gas  phase  reactions  to  start  playing  a  more  dominant  role. 
This  ratio  reaches  a  steady  value  of  about  10%  at  t  -  10s.  Heterogeneous  chemi¬ 
cal  kinetics  is,  thus,  indispensable  to  ignition,  feeding  back  into  gas  phase 
the  heat  necessary  to  raise  its  temperature  and  make  secondary  gas  phase  oxida¬ 


tion  possible 
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FIGURE  4. 


OVERALL  GAS  PHASE  TO  SURFACE  CONVERSION 
RATIO  DURING  THE  TRANSIENT 


Detailed  results 


Fig.  (5)  shews  the  axial  variation  of  the  ratio  of  fuel  conversion  by  cata¬ 
lytic  reaction  to  that  by  gas  phase  oxidation.  The  quantity  being  plotted  is: 


r  (x)  = 
gs 


fu  da  Ax 
aJ  g  _ 


„f  io  dl  Ax 
lJ  s 


(14) 


where  o  :  area  of  the  channel 

i  :  perimeter  of  the  channel 
Ax  :  the  local  cell  length. 

The  numerator  is  the  net  homogeneous  conversion  rate  at  an  axial  location 
while  the  denominator  is  the  local  heterogeneous  conversion  rate.  Initially, 
the  mixture  is  very  lean  and  conversion  of  fuel  is  restricted  to  the  surface  of 
the  catalyst.  This  trend  changes  as  the  incoming  fuel  mixture  turns  richer  and 
the  catalyst  heats  up.  The  contribution  of  homogeneous  kinetics  increases  more 
rapidly  near  the  inlet  of  the  channel  where  the  temperature  of  the  catalyst 
rises  more  rapidly  and  heat  is  transferred  to  the  gas.  As  the  temperature  of 
the  rear  end  comes  up,  gas  phase  contributions  increase  there,  too.  Finally,  at 
steady  state,  homogeneous  and  heterogeneous  kinetics  become  conparable  near  the 
midpoint  of  the  channel.  However,  very  little  conversion  takes  place  in  abso¬ 
lute  terms  at  this  location.  Most  of  the  fuel  is  used  up  near  the  inlet,  and 
where  the  homogenerous  reactions  still  play  a  minor  role. 

Figs.  (6a),  (6b)  and  (7)  provide  a  description  of  the  axial  redistribution 
of  thermal  energy  during  the  transient.  For  it,  the  following  quantities  are 
defined: 


4cataly tic 


+  +  Sp  +  Si 


(15) 


9.183E+00 
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where  9cataiytic  !  local  rate  of  catalytic  heat  release 

qg  :  local  rate  of  heat  storage  in  the  substrate 

q„  :  local  rate  of  heat  diffusion  into  the  gas  phase 

qR  :  local  rate  of  radiant  heat  loss 

qL  :  local  rate  of  conductive  radial  heat  loss  from 

the  substrate 

orQs  +  CG  +  CR  +  CL*51'  upon  normalizing  by  qcatalytic* 

In  Fig.  (6a)f  Qg  is  shown  at  different  times.  Initially,  the  profile  is  a 

straight  line  with  ~  1 •  This  implies  that  all  the  heat  released  goes  into 
the  substrate.  With  time,  the  axial  profiles  readjust.  While  near  the 
entrance,  the  profile  comes  down  to  Qs  ~  0,  toward  the  exit,  it  rises  so  that  ys 
>>  1.  This  indicates  that  the  front  end  of  the  catalyst  has  equilibrated  while 
the  rear  gets  heated  by  the  hotter  gas  (the  gas  gets  heated  by  the  solid  near 

the  entrance).  As  time  advances,  the  entire  axial  profile  is  seen  to  gradually 

return  toward  the  dotted  line,  representing  the  steady  state. 

In  Fig.  (6b),  i>G,  the  normalized  heat  flux  into  the  gas  is  shown  and  is 
seen  to  be  almost  a  mirror  image  of  Fig.  (6a).  Initially,  heat  transfer  into 
the  gas  is  negligible  and  Qq  ~  0.  A  period  of  readjustment  corresponding  to 
that  for  the  substrate  begins.  During  this  period,  Qq  ~  1  near  the  inlet  which 
indicates  that  most  of  the  energy  released  goes  towards  heating  the  gas.  Toward 
the  exit,  Qg  Is  a  large  negative  number.  This  is  a  direct  evidence  that  the  gas 
heated  by  the  solid  near  the  entrance,  travels  downstream  and  then  heats  the 
solid.  Finally,  the  profile  converges  toward  the  dotted  line,  representing  the 


steady  state 
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FIGURE  6a.  AXIAL  PROFILES  OF  HEAT  FLUX  FORM  THE  INTERFACE  TO  THE  SUBSTRATE 
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NONDIMENSIONAL  AXIAL  POSITION  (x/L) 


Radiative  losses  rise  rapidly  with  time  but  decrease  sharply  with  the 
distance  from  the  inlet,  as  seen  in  Fig.  (7).  Due  to  the  large  L/D  ratio  of  the 
channel,  radiation  ceases  to  be  of  significance  after  2  or  3  diameters  into  the 
channel,  and  only  the  first  10%  of  the  channel  is  shown.  A  similar  situation  is 
present  at  the  rear  end  of  the  catalyst.  However  the  associated  time  scale  is 
different  since  the  rear  end  of  the  catalyst  takes  a  long  time  heating  up. 

Also,  due  to  the  rear  end  of  the  test  section  getting  heated  by  the  exhaust 
gases,  the  fractional  radiation  exchange  is  significantly  lower  in  comparison  to 
that  at  the  inlet. 

Fig.  (8)  shows  axial  temperature  profiles  at  different  radial  locations, 
from  the  gas  centerline  to  the  substrate  half -thickness,  and  Fig.  9  shows  radial 
temperature  profiles  at  different  axial  locations.  Two  time  scales  can  be  iden¬ 
tified.  The  shorter  one  is  the  time  required  by  the  first  tenth  of  the 
substrate  to  reach  its  steady  state.  The  time  required  is  about  0.1s  which  is 
of  the  same  order  as  the  period  of  the  fuel  ramp.  It  suggests  that  near  the 
entrance  the  temperature  will  equilibrate  as  fast  as  the  fuel  is  supplied.  The 
other  time  scale  is  the  substrate  equilibration  time  and,  as  seen  earlier,  it  is 
about  10s.  The  characteristic  time  for  axial  heat  diffusion  in  the  substrate 
(Table  2)  is  of  the  order  of  103s.  The  significant  reduction  in  the  substrate 
heat  up  time  is  due  to  the  heating  of  the  back  part  of  the  solid  by  the  gas  and 
to  the  presence  of  axial  heat  sources. 

Comparison  with  steady  state  measurements 

No  transient  data  are  available  for  comparison  with  our  computations. 
However  the  steady  state  computed  through  the  transient  can  be  compared  with 
available  steady  state  data  [23].  It  is  realized  that  for  the  validity  of  the 
transient  model  satisfactory  agreement  with  steady  state  data  is  necessary  but 


not  sufficient 


NONDIMENSIONAL  AXIAL  POSITION  (x/L) 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

10 

o 

in 

o 

in 

o 

10 

o 

10 

o 

10 

(N 

T— 

CM 

CM 

T” 

T— 

(»)  3aruva3dW3i 


FIGURE  8. 


AXIAL  TEMPERATURE  PROFILES 
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FIGURE  9. 


RADIAL  TEMPERATURE  PROFILES 
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Comparisons  were  made  with  cases  differing  in  equivalence  ratio  and  mix¬ 
tures  velocity,  keeping  the  inlet  temperature  constant.  Computed  and  measured 
wall  temperatures  compare  favorably  as  shown  in  Fig.  (10),  and  so  do  predicted 
and  measured  exhaust  emissions  (Fig.  11).  A  detailed  discussion  of  the  steady 
state  can  be  found  in  Ref.  23.  Overall  the  results  indicate  the  predominance  of 
diffusion  rather  than  kinetics. 

CONCLUSIONS 

Following  are  the  main  conclusions. 

The  initial  heat  release  occurs  near  the  entrance  at  the  interface,  is 
controlled  by  heterogeneous  reactions,  and  goes  to  heat  the  solid  locally. 

The  temperature  of  the  solid  near  the  entrance  achieves  almost  its  final 
value  before  significant  heating  of  the  rest  of  the  solid  starts.  The  final 
temperature  of  the  solid  near  the  entrance  is  strongly  influenced  by  radiation 
losses.  Large  spatial  and  temporal  temperature  gradients  occur  in  the  solid  near 
the  entrance  controlled  mostly  by  the  availability  of  fuel.  These  gradients  may 
be  significant  as  far  as  fatigue  due  to  thermal  stresses  is  concerned. 

The  gas  is  heated  first  in  the  boundary  layer  near  the  entrance  where 
significant  radial  and  axial  gradients  are  present. 

The  gas  heated  near  the  entrance  flows  downstream  and  heats  the  back  of  the 
solid.  This  process  and  surface  reactions,  rather  than  conduction  along  the 
solid,  are  the  main  mechanisms  through  which  the  back  is  heated. 

Homogeneous  reactions  start  being  significant  first  in  the  boundary  layer 
near  the  entrance,  but  for  the  conditions  of  our  study  their  contribution  to  the 
overall  oxidation  of  CO  is  small  in  comparison  to  that  of  heterogeneous  reac¬ 
tions. 

Heterogeneous  oxidation  of  CO  soon  becomes  diffusion  controlled  even  though 
the  entire  light  off  transient  is  initially  controlled  by  heterogeneous  reaction 


U*  10  m/s 
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The  overall  equilibration  process  is  governed  by  the  thermal  inertia  of  the 
substrate  and  by  forced  convection.  99%  of  steady  state  is  achieved  in  about  10 
s  at  a  gas  velocity  of  10  m/s  and  in  1.5  s  at  70  m/s  for  the  particular 
substrate  substance  and  geometry  that  was  studied.  The  quasi-steady  gas  phase 
approximation  was  acceptable  for  our  cases  but  only  marginally  so  and  should  not 
be  assumed  to  be  valid  for  most  honeycomb  monoliths  of  practical  interest. 

Radiation  significantly  modifies  both  transient  and  steady  state,  par¬ 
ticularly  near  the  entrance  but  also  near  the  exit,  even  though  the  energy 
involved  is  small  and  because  of  its  localized  concentration. 


Perovskite  Catalysts  for  High  Temperature  Catalytic  Combustion* 

The  realization  of  catalytic  combustion  in  practical  combustion  systems  has 
in  part  been  limited  by  the  inadequacies  of  currently  available  combustion  cata¬ 
lysts.  One  of  the  objectives  of  Princeton's  Catalytic  Combustion  Program  has 
been  the  development  of  combustion  catalysts  with  improved  high  temperature  (> 
1500°K)  stability  in  terms  of  erosion,  corrosion,  and  thermal  stress  charac¬ 
teristics  and  improved  performance  characteristics  in  terms  of  lower  ignition 
temperatures,  extended  flammability  and  blowout  limits,  and  wider  fuel 
tolerances. 

In  this  section  results  are  presented  of  a  study  of  the  high  temperature 
durability  and  low  temperature  ignition  characteristics  of  platinum  doped 
perovskite  catalysts.  Operating  temperatures  of  catalytic  combustors  often 
exceed  1 500°K,  resulting  in  unacceptably  short  lifetimes  for  standard  catalysts. 
A  good  example  of  this  is  the  case  of  platinum  catalysts  which  have  very  good 
low  temperature  light  off  characteristics,  however  when  such  catalysts  are 
operated  at  typical  combustor  temperatures,  e.g.  1500°K,  the  platinum  coating  is 
lost  through  vaporization.  The  combustor  will  continue  to  operate  normally 
since  at  these  elevated  temperatures  the  ceramic  substrate  will  be  sufficiently 
catalytic  so  as  to  make  the  combustor's  performance  diffusion  controlled; 
however,  when  the  combustor  is  shut  off,  reignition  at  the  same  operating  con¬ 
ditions  will  not  be  possible  since  the  onset  of  ignition  is  controlled  by  the 
low  temperature  surface  reactions. 

A  new  catalyst  has  been  designed  to  exhibit  long  lifetime  at  high  tem¬ 
peratures  and  adequate  ignition  characteristics  at  low  temperatures.  The  cata¬ 
lyst  is  a  modified  perovskite  based  on  La (Cro. 5AI0. 5)^3.  This  ceramic  has  been 

*  Part  of  this  work  was  presented  at  the  Eastern  States  Section  of  the 
Combustion  Institute  Fall  Meeting,  December  1982. 
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doped  to  make  it  electrically  conductive,  and  consequently  it  can  be  resistively 
heated  to  bring  the  catalyst  up  to  the  required  light-off  temperature.  In  addi¬ 
tion,  platinum  has  been  incorporated  into  the  crystal  structure  to  give  improved 
low  temperature  light-off  while  having  a  low  platinum  vapor  pressure  at  high 
operating  temperatures. 

This  catalyst  can  be  used  in  powdered  form,  by  washcoating  it  onto  a  high 
temperature  ceramic  substrate,  or  as  sintered  monolithic  structures,  e.g. 
plates.  The  advantage  of  using  the  catalyst  in  the  form  of  plates  is  that  they 
can  be  resistively  heated  to  assist  light-off.  However,  the  catalyst  is  more 
readily  available  in  powdered  form  and  therefore  the  first  tests  with  the  new 
catalyst  have  been  made  with  it  washcoated  on  a  honeycomb  substrate.  The 
substrate  used  was  mullite,  three  inches  long  with  1/16  inch  square  cells. 

Experiments  using  a  platinum  washcoat  and  perovskite  powder  washcoats  with 
nominal  0.1%  (by  weight)  platinum  (supplied  by  Transtech,  Inc.)  and  with  nominal 
0.5%,  1.0%  and  2.0%  (by  weight)  platinum  (supplied  by  Prof.  Harlan  Anderson, 
University  of  Missouri -Rolla )  have  been  conducted.  The  tests  consisted  of 
establishing  the  light-off  temperature  and  low  temperature  performance  for 
several  equivalence  ratios  and  inlet  velocities,  after  which  the  catalyst  was 
"aged"  for  several  hours  at  1500°K  and  then  the  light-off  and  low  temperature 
experiments  repeated.  The  substrate  temprature  profile  for  a  typcial  aging  run 
is  shown  in  Figure  12.  The  purpose  of  running  at  low  temperature  is  to  insure 
that  all  of  the  heat  release  is  due  to  surface  reactions.  Therefore  the  maximum 
substrate  temperature  in  these  runs  was  kept  below  1050°K.  Even  at  low  tem¬ 
peratures  however,  it  is  important  that  the  overall  process  be  surface  reaction 
rate  controlled  and  not  diffusion  controlled.  Then  the  measured  substrate  tem¬ 
perature  profiles  will  provide  an  indication  of  changes  in  surface  activity. 

Figure  13  shows  the  substrate  temperature  profiles  before  and  after  aging 
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FIGURE  13.  PURE  PLATINUM  CATALYST  PROPANE  FUEL,  QUIVALENCE  RATIO 
INLET  TEMP.  -  589K,  INLET  VELOCITY  *  10.1  m/s 


for  the  platinum  coated  catalyst  using  propane  fuel  at  an  equilvalence  ratio  of 
0.28,  inlet  temperature  of  580°K,  inlet  velocity  at  10  m/sec  and  1  atmosphere 
pressure.  After  aging  for  four  hours  at  1500°K  the  catalyst  would  not  reignite 
at  these  operating  conditions,  clearly  indicating  a  loss  of  surface  activity 
which  in  this  case  is  due  to  a  loss  of  platinum  as  a  result  of  vaporization. 

The  results  for  the  same  catalyst  before  and  after  aging  using  hydrogen  fuel  are 
shown  in  Figure  14.  In  this  case  for  the  same  catalyst  there  is  no  indication 
of  a  change  in  surface  activity.  The  reason  for  this  is  that  hydrogen  is  so 
reactive  on  platinum  that  the  catalyst's  performance  is  diffusion  controlled 
even  after  substantial  platinum  loss.  These  two  cases  are  interesting  in  that 
the  first  clearly  demonstrates  the  problem  of  catalyst  aging  and  the  second 
demonstrates  the  importance  of  evaluating  the  catalyst's  surface  activity  under 
surface  reaction  rate  controlled  conditions. 

Figure  15  shows  measured  substrate  temperature  profiles  for  the  perovskite 
with  nominal  1%  by  weight  platinum  using  propane  fuel  at  an  equivalence  ratio 
of  0.3,  inlet  temperature  of  700°K,  inlet  velocity  of  6  m/sec,  and  1  atmosphere 
pressure.  The  before  aging  result  is  with  the  new  catalyst.  Also  shown  are 
results  after  four  successive  aging  runs.  The  indicated  times  are  cumulative. 

As  the  results  show,  there  is  a  noticeable  loss  of  activity  after  the  second  and 
fourth  aging  runs  which  were  done  with  a  mixture  of  hydrogen  and  propane. 
However,  after  the  first  and  third  aging  runs,  which  were  done  with  pure  pro¬ 
pane,  there  is  little  indication  of  reduced  surface  activity.  These  results  are 
very  encouraging,  indicating  very  little  change  in  surface  activity  when  used 
with  propane  fuel.  The  tests  with  hydrogen  however  indicate  a  definite  reduc¬ 
tion  in  the  activity  of  the  catalyst. 

In  order  to  confirm  these  results  plans  were  made  to  repeat  the  tests  with 
the  platinum  doped  perovskite  using  only  propane  during  the  aging  runs.  Three 
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FIGURE  14.  PURE  PLATINUM  CATALYST  HYDROGEN  FUEL,  EQUIVALENCE  RATIO 
INLET  TEMP.  -  590  K,  INLET  VELOCITY  -  10.2  m/s 
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FIGURE  15.  PEROVSKITE  WITH  NOMINAL  1%  BY  WEIGHT  PLATINUM  CATALYST 
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new  catalysts  were  prepared  with  new  perovskite  powders  (supplied  by  Prof.  H. 
Anderson)  having  0.5%,  1.0%  and  2.0%  by  weight  platinum  content.  A  perovskite 
powder  with  0.1%  by  weight  platinum  supplied  previously  by  Transtech  Inc.  was 
also  tested.  The  measured  substrate  temperture  profiles  with  these  four  cata¬ 
lysts  are  shown  in  Figure  16  using  propane  fuel  at  an  equivalence  ratio  of  0.3, 
inlet  temperature  of  800°K,  inlet  velocity  of  10  m/sec,  and  1  atmosphere 
pressure.  Also  shown  for  comparison  is  the  corresponding  temperature  profile 
with  the  original  1%  platinum  catalyst.  This  comparison  clearly  shows  that  all 
four  of  the  new  catalyst  samples  were  considerably  less  active  than  the  original 
catalyst  even  though  in  one  case  the  indicated  platinum  content  was  twice  as 
great. 

The  remainder  of  this  report  is  a  summary  of  the  status  of  the  investiga¬ 
tion  of  the  difference  in  surface  activity  between  the  new  perovskite  powders 
and  the  original,  considerably  more  active  1%  platinum  doped  perovskite  powder. 
The  first  test  was  basically  a  confirmation  of  the  flow  reactor  results.  This 
involved  testing  the  perovskite  catalysts  using  a  differential  scanning  calori¬ 
meter  (DSC).  This  simply  involves  heating  two  small  samples,  one  of  the  cata¬ 
lyst  and  the  other  an  inert,  in  a  temperature  programmed  oven  and  monitoring  the 
temperature  difference  between  the  catalytic  and  the  inert  samples  which  is  an 
indication  of  the  heat  release  due  to  catalytic  reactions.  The  oven  can  be 
filled  with  any  fuel/air  mixture  to  test  the  catalytic  activity  with  respect  to 
fuel  identity  and  equivalence  ratio.  In  addition  to  the  original  1%  sample  and 
the  new  0.1%,  0.5%,  1.0%  and  2.0%  samples,  tests  were  run  with  platinum  black 
and  a  perovskite  powder  with  no  platinum.  The  results  from  the  DSC  tests,  run 
with  propane  in  air  at  an  equivalence  ratio  of  2.8,  are  shown  in  Figure  17  and 
basically  confirm  the  flow  reactor  results.  The  temperature  at  which  the  heat 
release  drastically  increases  can  be  thought  of  as  a  relative  ignition  tern- 
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perature.  Note  that  the  relative  ignition  temperatures  of  the  new  catalysts  are 
consistent  with  their  relative  platinum  content  but  that  the  original  1%  sample 
is  considerably  more  active,  i.e.  has  a  lower  ignition  temperature.  In  fact  the 
original  1%  catalyst  is  nearly  as  active  as  the  platinum  catalyst.  An  important 
outcome  of  the  DSC  testing  was  the  demonstration  that  such  a  method  could  be 
used  to  determine  the  relative  activity  of  catalyst  samples  with  various 
fuel/air  combinations.  This  has  proven  to  be  an  efficient  and  convenient  method 
of  screening  new  catalysts  before  preparing  a  washcoated  substrate  and  testing 
in  the  flow  reactor. 

After  the  differences  in  surface  activity  were  confirmed  with  the  DSC, 
several  tests  were  made  to  establish  the  cause  of  this  difference.  The  first  of 
these  tests  involved  using  X-ray  fluorescence  to  determine  the  actual  platinum, 
chronium  and  lanthanum  content  of  the  perovskite  catalyst  samples.  The  results 
from  these  measurements  are  given  in  Tables  4,  5  and  6  where  the  absolute  plati¬ 
num  composition  and  the  relative  chromium  and  lanthanum  compositions  are  given. 
The  measured  platinum  content  was  found  to  agree  reasonably  well  with  the 
ejected  platinum  content.  The  small  differences  cannot  explain  the  drastic 
differences  in  catalyst  activity.  The  relative  chromium  and  lanthanum  com¬ 
positions  were  found  to  be  nearly  the  same  in  all  samples  and  again  the  small 
differences  cannot  explain  the  differences  in  catalyst  activity. 

The  second  test  involved  scanning  electron  microscope  studies  of  the 
crystal  morphology  and  sintering.  Photographs  of  the  original  1%  and  the  new 
0.5%,  1.0%  and  2.0%  catalysts  are  shown  in  Figure  18  at  a  magnification  of  5400. 
The  basic  conclusion  from  these  photographs  is  that  there  is  no  difference  in 
particle  size,  shape  or  degree  of  sintering  that  can  explain  the  observed  dif¬ 
ference  in  catalytic  activity. 


The  third  test  was  BET  measurements  of  the  relative  surface  area.  The 
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TABLE  4  -  PLATINUM  CONCENTRATION 


SAMPLE 

MEASURED 

• 

CONCENTRATION 

0.1%  Pt 

(8/81  ) 

0.1% 

0.5%  Pt 

(11/82) 

0.34% 

• 

1.0%  Pt 

(11/82) 

1.03% 

2.0%  Pt 

(11/82) 

2.32% 

1.0%  Pt 

(8/82) 

0.98% 

TABLE  5  - 

■  RELATIVE 

LANTHANUM  CONCENTRATION 

tf 

SAMPLE 

MEASURED  AREA 

0.5%  Pt 

(11/82) 

9816 

< 

1.0%  Pt 

(11/82) 

9677 

2.0%  Pt 

(11/82) 

9867 

1.0%  Pt 

(8/82) 

9861 

TABLE  6  - 

RELATIVE 

CHROMIUM  CONCENTRATION 

• 

SAMPLE 

MEASURED  AREA 

0.5%  Pt 

(11/82) 

45693 

1.0%  Pt 

(11/82) 

48315 

Q 

2.0%  Pt 

(11/82) 

49376 

1.0%  Pt 

(8/82) 

45671 

TABLE 

7  -  BET  RELATIVE  SURFACE  AREA 

SAMPLE 

SURFACE  AREA 

0.5%  Pt 

(11/82) 

8.073  m2/g 

\ 

1.0%  Pt 

(11/82) 

10.115  m2/g 

2.0%  Pt 

(11/82) 

6. 876  m2/g 

1.0%  Pt 

(8/82) 

9.552  m2/g 

< 


results  of  these  tests  are  given  in  Table  7.  Although  there  are  differences  in 
the  effective  surface  area  they  do  not  correlate  with  the  surface  activity 
differences. 

The  fourth  and  final  test  to  date  involved  the  use  of  Fourier  transform 
infrared  spectroscopy  (FTIR)  to  characterize  the  catalyst.  FTIR  spectra  of  the 
original  1%  catalyst  and  the  new  2%  catalyst  are  shown  in  Figure  19.  Although 
the  spectral  characteristics  have  not  yet  been  identified  there  are  definite 
differences  in  the  two  spectra,  particularly  near  3400  cm"'  and  1500  cm-',  which 
may  be  related  to  the  observed  variations  in  catalytic  activity.  Work  is  con¬ 
tinuing  to  identify  the  spectral  characteristics  which  in  turn  may  explain  the 
observed  differences  in  catalytic  activity. 

In  summary,  perovskite  catalysts  with  small  amounts  of  platinum  have  been 
shown  capable  of  achieving  low  temperature  light-off  and  high  temperature  stabi¬ 
lity  in  preliminary  tests.  Attempts  to  reproduce  these  results  have  revealed 
unexplainable  variations  in  catalytic  activity  from  sample  to  sample. 
Measurements  of  platinum,  lanthanum  and  chromium  composition,  of  crystal  morpho¬ 
logy  and  sintering,  and  or  relative  surface  area  have  not  provided  an  explana¬ 
tion  for  these  differences  in  catalytic  activity.  Preliminary  infrared  spectra 
of  the  catalyst  do  show  differences  which  may  or  may  not  be  related  to  the 
observed  variations  in  catalytic  activity.  Work  is  continuing  to  identify  the 
spectral  characteristics  of  the  different  perovskite  samples,  which  in  turn  may 
provide  insights  to  explain  the  differences  in  catalytic  activity. 

The  idea  of  evaluating  the  relative  performance  of  the  catalysts  in  simple 
bench  top  experiments  before  actually  testing  them  in  the  more  complicated  flow 
reactor  experiments  is  very  attractive.  As  mentioned  earlier,  the  differential 
scanning  calorimeter  (DSC)  has  proven  to  be  a  simple  and  convenient  method  of 
characterizing  the  relative  catalytic  activity.  In  the  final  study  covered  in 
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this  report  the  DSC  was  used  to  characterize  the  "aging"  of  the  2%  Pt-perovskite 
powder  before  and  after  exposure  to  high  temperature  in  both  oxygen  and  nitrogen 
environments.  The  samples  were  heated  in  a  tube  furnace  to  temperatures  as  high 
as  1 500°K  for  periods  of  4  hours.  In  Figure  20  the  DSC  plots  for  the  new  and 
various  aged  samples  are  given.  As  3hown  the  catalytic  activity  drastically 
deteriorates  in  the  high  temperature  oxygen  environment  whereas  comparable  aging 
in  nitrogen  shows  no  appreciable  change  in  catalytic  activity.  Also  shown  is 
that  the  oxygen  exposed  catalysts  can  be  substantially  regenerated  in  a  high 
temperature  oxygen-free  environment.  At  this  point  the  effect  of  the  oxygen  on 
the  catalyst  is  not  understood. 

Fourier  transform  infrared  spectra  of  the  new  and  variously  aged  2%  Pt- 
perovskite  catalysts  were  also  taken  and  are  shown  in  Figure  21.  There  is 
however  no  apparent  correlation  between  the  infrared  spectra  and  the  observed 
catalyst  aging  phenomena. 

At  the  conclusion  of  the  grant  period  covered  by  this  report  work  was  con¬ 
tinuing  to  determine  the  cause  of  these  variations  in  catalytic  activity. 
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surface  area  of  the  catalyst,  m2 
cross-sectional  area  of  the  substrate,  m2 
radial  heat  loss  constant 

specific  heat  at  constant  pressure,  J/kg-K 

binary  diffusion  coefficient,  m2/s 

hydraulic  diameter  of  the  channel, m 

activation  energy  for  homogeneous  and  heterogeneous 

reactions,  kcal/mole 

enthalpy  per  unit  mass,  J/kg 

pre-exponent  for  the  homogeneous  reaction  rate, 
cm/mole-s 

pre-exponent  for  the  heterogeneous  reaction  rate, 
cm/s 

thermal  conductivity,  J/m-S-K 
perimeter  of  the  channel,  m 
length  of  the  channel,  m 

CO  mass  consumption  rate  per  unit  area,  kg/m2-s 

Nusselt  number 

pressure,  N/m2 
heat  flux,  W/m2 

heat  release  rate  due  to  heterogeneous  reaction  per 
unit  area,  W/m2 

normalized  heat  transfer  rate,  or 

heat  of  combustion  of  homogeneous  reaction,  J/kg 

channel  radial  coordinate,  m 

channel  radius,  m 


universal  gas  constant  cal/mole-K 
thickness  of  the  substrate,  m 
source  term 
time,  s 

temperature,  K 
axial  velocity,  m/s 
inlet  velocity,  m/s 
radial  velocity,  m/s 
volume  of  channel  gas,  m3 
molecular  weight  of  species  k 
channel  axial  coordinate,  m 
mole  fraction  of  species  k 
mass  fraction  of  species  k 

thermal  diffusivity,  m2/s,  or  absorptivity 

transport  property 

emissivity 

gas  viscosity,  kg/m-s 
stoichiometric  coefficient 
density,  kg/m3 

cross-sectional  area  of  the  channel,  m2»  or  stress 
component,  N/m 2,  or  Stefan  Boltzmann  constant, 
W/m2-*4 

non-dimensional  time  or  characteristic  time,  s 
any  dependent  variable  (such  as  u, v, p, . . .etc. ) ,  or 
equivalence  ratio 


reaction  rate 
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in, IN, 0 


out, OUT, L 


Superscript 


backward  reaction 


forward  reaction,  or  face  of  the  substrate,  or  final 


inlet,  i.e.  x  ■  0 


species  index 

1-th  homogeneous  reaction 


outlet,  i.e.  x  =  L 


substrate 


test  section 


wall  (substrate) 


reactant 


product 


substrate  side  of  the  interface 


gas  side  of  the  interface 
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